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Abstract 

The single photon annihilation contributions for the positronium ground state hyperfine 
H ' splitting are calculated analytically to order mecfi using NRQED. Based on intuitive physical 

arguments the same result can also be determined by a trivial calculation using results from 
existing literature. Our result completes the hyperfine splitting calculation to order meOt^ ■ 
We compare the theoretical prediction with the most recent experimental measurement. 



Positronium, a two-body bound state consisting of an electron and a positron, belongs to the first 
systems studied within the relativistic quantum theory developed by Dirac. The existence of positro- 
nium was predicted in 1934 and experimentally verified at the beginning of the 1950s [^. For the 
ground state hyperfine splitting (hfs), the mass difference between the 1^5i and I^Sq state, steadily 
improved experimental measurements have meanwhile reached a precision of 3.6 ppm which makes 
the calculation of all 0(0^) corrections to the leading and next-to-leading order expression mandatory. 
So far only the order a^lna"^ corrections have been fully determined Q. Including also the known 
0{a'^ In^ a~^) corrections ||5|, |^ the theoretical expression for the hfs readsQ 
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We have calculated the single photon annihilation (1 — 7 ann) contributions to the constant K. All 
other contributions to K coming from the non-annihilation, and the two and three photon annihilation 
processes have been calculated before (see Table |l]). Although our result completes the calculations to 
order m(,a^ the theoretical situation remains unresolved due to a discrepancy in two older calculations 
for some of the non-annihilation contributions to order 171^0^ ■ 

In this letter we report on two methods to determine Wnnlo""- The first one is systematic using 
NRQED 1^ , which is based on the concept of effective field theories and the separation of effects from 
non-relativistic and relativistic momenta, and the second one relies on physical intuition using two 
results from already existing literature. We also summarize the status of the theoretical calculation 
to the hfs in view of the most recent available experimental data. We would like to note that the 
presentation of our NRQED calculation is only meant to illustrate the basic steps of our calculation. 
A more detailed and explicit work on the NRQED method is in preparation. Our second method, on 
the other hand, is almost trivial and represents a true "back of the envelope" calculation. 

For the NRQED calculation we start from the NRQED Lagrangian Q 
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where if) and x ^^'s the electron and positron Pauli spinors and Dt and D the time and space compo- 
nents of the gauge covariant derivative -D^. In Eq. (|2|) the straightforward bilinear positron terms are 
omitted and only those four fermion interaction relevant for the one photon annihilation contributions 
to the hfs are displayed. The renormalization constants ci, . . . , C4, di, ^2 are normalized to one at the 
Born level. For Wnnlo"" only the radiative corrections to di have to be calculated. To order m-eO^, and 
if we consider only the contributions from the one photon annihilation graphs, all retardation effects 
can be neglected. This means that the transverse photon propagators can be used in the instanta- 
neous approximation, i.e. their energy dependence is dropped. Indeed, simple counting rules show 
that retardation corrections to the one photon annihilation diagrams would set in at order meoJ . In 
the instantaneous approximation, all NRQED interactions can be written as two-body potentials. In 



^ We use natural units, in which fi — c = 1. 
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Figure 1: Interaction potentials contributing to VFj^nlo'"- Kin denotes the relativistic kinetic energy 
correction coming from the D'^/Sm^ terms in the NRQED Lagrangian. 



momentum space representation all the potentials needed for the present calculation are given by (see 
Fig. 0) 
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where A is a small fictitious photon mass to regularize IR divergences and Szf denotes the elec- 
tron/positron spin operator. V^^i denotes the l/m-g corrections to the Coulomb potential including 
longitudinal and transverse photon exchange. V4 accounts for the leading order annihilation process 
e+e" ^ 7 — > e+e^ and V4der denotes relativistic corrections to V4 from the energy dependence of the 
annihilation photon and from the l/mg contributions in the Dirac spinors. The calculation of 
proceeds in two basic steps. 
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1. Matching calculation: calculation of the 0{a) and 0(02) contributions to the constant di by 
matching the NRQED and QED amplitudes for the elastic s-channel scattering of free and on- 
shell electrons and positrons via a single photon, e+e~ — > 7 ^ e+e", up to two loops and to 
NNLO in the velocity of the electrons and positrons in the cm. frame. 

2. Bound state calculation: calculation of W^'^j^q'' by solving the non-relativistic bound state prob- 
lem in form of the Schrodinger equation (i.e. including the non-relativistic kinetic energy and the 
Coulomb interaction) exactly and by treating the relativistic effects using Rayleigh-Schrodinger 
time- independent perturbation theory (TIPT). 

Matching calculation: to determine the NRQED amplitude e^e~ — > 7 — > e^e~ for free and on-shell 
electrons and positrons up to two loops and NNLO in the velocity the diagrams displayed in Fig. 
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Figure 2: NRQED diagrams for the matching calculation. Permutations are not displayed. 



have to be calculated. It is sufficient to consider only scattering of the e~^e~ pair in a ^Si state 
because a ^Sq state cannot annihilate into a single photon due to C invariance. We have regularized 
all UV divergent integrations by using a momentum cut-off. As a consequence the finite terms in 
the NRQED amplitude (and also in the constant di) depend on the routing of the loop momenta 
through the diagrams |^]. Therefore, to be consistent exactly the same way of routing has to be 
used in the bound state calculation. We come back to this point later. The corresponding QED 
amplitude for the scattering process has to be determined by using conventional covariant multi-loop 
perturbation theory. Whereas the one-loop results for the vertex corrections [^] and the one and 
two-loop contributions to the vacuum polarization function |llO| have been known for quite a long 



time, the two-loop vertex corrections have been calculated recently by one of the authors |11|. The 
QED amplitude is renormalized by the common on-shell renormalization scheme, where a is the fine 
structure constant and the wave function renormalization constant is fixed by the requirement that the 
residue of the fermion propagator is one. The 0{a) and O(a^) contributions of the renormalization 
constant di are then determined by demanding equality of the NRQED and QED amplitude at the 
one- and the two-loop level. Because all IR divergent and velocity dependent contributions are equal 
in NRQED and QED, di contains only UV divergent and constant contributions. 
Bound state calculation: to finally determine VI^nnlo'" we start from the well known solution of the 
non-relativistic positronium problem (in form of the Schrodinger equation) and incorporate V4, V4dcr) 
V^ci and via first and second order TIPT. For each insertion of V4, the contributions from di have 
also to be taken into account. The divergences in di automatically remove the UV divergences which 
arise in the bound state calculation. At this point we want to emphasize again that, to be consistent, 
the routing of the momenta in the bound state calculation has to be exactly the same as the routing 
in the NRQED scattering diagrams. Also the finite terms in the bound state integrals depend on the 
routing. Combining the result of the bound state integrals with the contributions in di leads to the 
cancellation of the routing-dependent terms. We have checked this fact by choosing different routings 
in our calculation. 
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The final result for the order rriga^ one photon annihilation contributions to the hfs readsQ 
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The lna~^ term was already known and is included in the Ina^^ contribution quoted in Eq. (||). The 
one photon annihilation contribution to the constant K corresponds to a contribution of —2.34 Mhz 
to the theoretical prediction of the hfs (see Table |l|) . 

The second, more intuitive method to determine Wj^NLo" starts from the formal result for 
the energy shift due to one photon annihilation for S-wave triplet bound states with radial quantum 
numbers n using first to third order TIPT^, 
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where \ i = l,m,k, represent normalized (bound state and continuum) eigenfunctions to the positro- 
nium Schrodinger equation with the eigenvalues Ei. It is evident from the form of the operator V4 that 
W^"^ ™ depends only on the zero-distance Coulomb Green function An = (0 | ^i^n e -_e I ^) (where 



the n'^Si bound state pole is subtracted) and on the rate for annihilation of a n^Si bound state into a 
single photon, Pn = {n\V4\n) + [{n\Vi J;^„ eI-Ei + ^^in) I ) + H.c] + ( n | Kdcr I n ) (where the 
effects from Kei and Kin are included in form of corrections to the wave function). Because An and 
Pn are UV divergent from the integration over the high energy modes, they have to be renormalized. 
In the NRQED approach this was achieved by the renormalization constant di. Here, renormalization 
will be carried out by relating An and Pn to physical (and finite) quantities which incorporate the 
proper short-distance physics from the one photon annihilation process. For An this physical quantity 
is just the QED vacuum polarization function in the non-relativistic limit and for P„ the abelian 
contribution of the NNLO expression for the leptonic decay width of a superheavy quark-antiquark 
n^Si bound state |13]. Both quantities have been determined recently in | 
of [14, |l^ it is straightforward to derive the renormalized versions of An and Pn 
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^ The contributions in Wn'nlo'" coming from the vacuum polarization effects of the annihilation photon have been 
calculated before in [ p^ The result in |l2j contains an error in the treatment of the one-loop vacuum polarization 
(see U4]). The vacuum polarization contributions calculated in jl^] are in agreement with our result. 

5*0 states do not contribute to VK^"^ because they cannot annihilate into a single photon. States with higher 
orbital angular momentum, on the other hand, are irrelevant because their wave functions vanish for zero electron-positron 
distance. 
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where Ci = 2^(-3 + ^Cs) - X5 + i In 2 and C2 = :^(^ - Ca) + | In 2 - f|. 7 is the Euler constant 
and ^ the digamma function. Inserting A'^^'' and P^^^" back into expression @ we arrive at 
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It is an easy task to check for the ground state n = 1 that in Eq. ( [ll| ) the well known order meCt^ 
[Wlo = J meO^) and mga^ (W^n'Jo''"" = — meO^) one photon annihilation contributions to the hfs 
are correctly reproduced and that the mea^ contributions are equal to TVnnlo'") Eq. (|^. This second, 
very simple intuitive method to determine TVnnlcd" does not only represent a nice cross check for the 
systematic NRQED calculation but also illustrates that W'^~~' is directly related to other physical 
quantities. Beyond order mga^ expression ([III) is not valid because essential retardation effects are 
not taken into account. 
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204 386.7(1) 
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19.1 
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Sum (Caswell-Lepage) 
Sum (Pachucki) 
Experiment 



203 388.3(6) 
203 392.2(1) 
203 389.1(7) 



Table 1: Summary of the theoretical calculations to the hfs. Only the references with the 
first correct calculations are given. 



In Table | we have summarized the status of the theoretical calculation to the hfs of the 
positronium ground state including our own result. To order meOt^ the logarithmic in a and constant 
contributions are given separately. The constant terms are further subdivided into non-annihilation, 
and one, two and three photon annihilation contributions. The error in the order meO^ result (1.) 
comes from the uncertainties in the input parameters a, fi and nif. and the errors in 4. and 5. 
are numerical. For all other contributions the errors are negligible. As indicated, there are two 
contradictory calculations for some of non- annihilation contributions based on results from Caswell 
and Lepage (C/L) Q and Pachucki (Pa) |20|. The result containing the Caswell-Lepage calculation 
leads to perfect agreement between theory and experiment (VFth — W^^ = —0.8(1.0) Mhz), whereas 
the hfs prediction based on the result by Pachucki leads to a discrepancy of more than four standard 
deviations (Wth — W^cx = 3.1(0.7) Mhz). It remains the task of future examinations to finally resolve 
the theoretical situation. 
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During completion of this work we were informed of work on the same subject by Adkins, Feh 
and Mitrikov using the Bethe-Salpeter formalism and numerical methods. Their result agrees with 
ours representing an independent cross check. We thank G. Adkins and his group for reporting their 
result to us prior to publication. We also thank G.P. Lepage for useful discussions. This work is 
supported in part by the Department of Energy under contract DOE DE-FG03-90ER40546 and by 
the Natural Sciences and Engineering Research Council of Canada. 
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